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I. INTRODUCTION 

Scalar-Tensor (ST) gravity theories constitute one of the main sources of extensions to General 
Relativity (GR) (for reviews see [DU]). They constitute a class of theories with well posed actions 
where the gravitational fields are given by both the metric and a scalar field. Diverse physical 
motivations for GR extensions naturally lead to ST actions m- On the other hand, currently 
there is no solid observational or experimental evidence that clearly favours a ST model over GR. 
This has lead to considerations on mechanisms that could somehow hide the non-observed scalar 
field effects in the Solar system, while allowing it to develop a nontrivial role at different scales, 
say at the cosmological level. Such mechanisms, which change the effective number of degrees 
of freedom of gravity, by hiding the scalar field at appropriate scales, have received considerable 
attention recently [Mi- 
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External fields appear in diverse contexts in physics. In particular, an external field is sometimes 
defined as a field whose expression as a spacetime function is known beforehand. Thence the role 
of the field equations is to provide solutions to the other (non-external) fields. From the action 
perspective, an external field is a quantity that should not be varied, since it is already known. 
For instance, the Minkowski metric in the context of electromagnetism in flat spacetime is an 
external rank two field. In this case the spacetime geometry is fixed beforehand and one uses the 
field equations, which are derived from variations with respect to other fields, in order to derive 
electromagnetic configurations that are compatible with such geometry. 

It is not a novelty that external fields are used as part of the framework of gravity theories. For 
instance, in GM3] the authors, based on Renormalization Group expectations, argue in favor of a 
gravitational action with two external scalar fields (G and A). This is a natural starting point since 
the dynamics of these two fields, in this context, should not be inferred from pure classical field 
equations, the Renormalization Group equations together with certain ansatz are used to partially 
fix the dynamics of G and A (for recent developments see jT2j) • Another proposal that uses external 
fields is that of massive gravity in arbitrary backgrounds. The latter is equivalent to bigravity in 
the limit in which the coupling constant of one of the rank-2 tensors goes to infinity (see e.g. m 
see also D3H32D- These actions contain a usual rank-2 tensor together with an additional one that 
has external character, that is the full field equations do not consider an equation coming from the 
action variation with respect to the external tensor (which may or may not describe a Minkowski 
spacetime geometry). We also cite the case of unimodular gravity (see [T8H23] and references 
therein). In unimodular gravity, when formulated as in [[183, there is no field that is fully external, 
but the determinant of the metric behaves as an external field. That is, the metric has the same 
role as in the Einstein-Hilbert action, except that its determinant has a fixed valuej^] Its action 
can also be recast by using a Lagrange multiplier such that the metric is fully non-external, but a 
condition is imposed between the metric determinant and an external scalar density field mm- 
As a last example, the more recently proposed approach to modify gravity inspired on continuous 
mechanics [24] , is closely related to the approach developed here. In Ref. [HJ an external tensor 
(which is not a metric) is also used. Contrary to Ref. [24] . our emphasis here will be on an external 
scalar and the comparison between the frameworks of GR and ST gravity to EST gravity, that is, 
a ST gravity whose scalar is an external field. 

In order to be clear, the definition that is here used for an external field is the following: it is 

1 In the context of unimodular gravity, it is common to say that the metric determinant is “non-dynamical” (e.g., 
mi), but the label of “external”, as used here, also applies to it and it has been used in some references, in 
particular see ES3- 
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a field whose variation inside the action should not be considered in the process of deriving the 
field equations. Hence, for a fixed number of fields, the higher is the number of external fields, the 
lower is the number of field equations derived from the action. The main concern in this work is 
precisely to look for EST solutions that are different from both GR and ST solutions. It turns out 
that EST theories allow a different way of recovering GR behaviour in some spacetime regions. 
The transition between GR to non-GR behaviour depends on the boundary conditions associated 
with the system. 

This work is organised as follows: Section [IT] is devoted both to the introduction of EST gravity 
and to address general differences and similarities to ST gravity and GR. A particular cosmological 
solution that exhibits both a GR phase and a ST one (but that is globally neither a GR nor a ST 
solution) is proposed in Section III In Section IV we present our conclusions and in Appendices [A| 
10 we present respectively issues on the stress-energy conservation related with external fields, 
additional details on the proposed particular solution and an extension of the latter for barotropic 
fluids. A mechanical example on the use of external variables is presented in Appendix [Dj and 
Appendix [E] shows a first test using SN data. 

Commonly used abbreviations: BD = Brans-Dicke, EST = Scalar-Tensor with an external 
scalar field, GR = General Relativity, RW = Robertson-Walker, ST = Scalar-Tensor. 


II. SCALAR-TENSOR GRAVITY WITH AN EXTERNAL SCALAR FIELD 

Consider the following action in which 4> is an external scalar field 

S[g,y] = ^-J (*R - ^l g - b d a <S>d b <S> - 24>A ($)) J-gd A x + S m [g, T], (1) 

where \k collectively represents matter fields. It looks as the usual ST action, but the difference (at 
the action level) is that the scalar field $ is not to be varied, since it is external (as defined in the 
Introduction). For this reason, the action dependence in <f> does not explicitly enter the notation 
l S'[g,\k]. Thence, apart from the matter T field equations, the only equations derived from the 
action are given by 

$ ( R ab - l -Rg ab + A g a ^j = 8vr T ab + | (d a ^ - \g ab g cd d c ^d d ^j + V a 9 b <h - g ab ,□$. (2) 

Since (i) 4> does not appear in the matter sector of the action (jT]) , (ii) all the T-dependent terms are 
in S m , and (iii) none of the fields entering S m is external, one gets the stress tensor conservation 


V a T ab = 0 


( 3 ) 
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as a direct consequence of the diffeomorphism invariance of S m and of the fact that each field 'I' 
generates a Lagrange equation (see Appendix [A] for the requirement of each assumption). 

Since <f> is external, there is no field equation derived from the action variation with respect 
to Thence, equations ([ 2 ]) and ([ 3 ]) are the only equations of the theory. Thanks to ([ 3 ]), the 
contracted Bianchi identity and the covariant derivatives commutation rule, the divergence of ([ 2 ]) 
leads to 


2^DT + + R ~ 2 7T ($ A ) 

$ d§ Vd)/ d<f> 

Contracting ([ 2 ]) and eliminating R, one gets 

87 tT - Jd c $d c <$> - 2d> (A - d*A') 


d b $ = 0 . 


- 


+ 3 


&<!> = 0 . 


(4) 


(5) 


We remark that this equation is not the scalar field equation that is valid in usual ST gravity]^] 
because of the presence of as an overall factor. In the case that <9<h 7 ^ 0 everywhere, then eq. 
© leads to the usual ST equation 

8 vrT - Jd c $d c <5> - 2<3? (A - <f>A') 


□<f> = 


2w + 3 


( 6 ) 


The above equation is precisely the same equation that is derived from the action ([Tj) variation 
with respect to T. It then should be clear that the equation set ©, © and ([6]) composes the 
usual ST fundamental field equations 12 EE while the equation set © and Q composes the EST 
fundamental field equations (from which ([ 5 ]), that is less constraining than ([h]), is just a logical 
consequence). Let us also stress that, in the usual ST case, the contracted Bianchi identity just 
results in a triviality, since equation ([b]) is already known as a fundamental held equation. On the 
other hand, in the EST case, the contracted Bianchi identity leads to equation ([ 5 ]), that is worth 
to be considered, since it does not appear explicitly (while it is implicitly present) in the set of 
fundamental held equations. Therefore, the EST action ©, together with the assumption that 
9T 7 ^ 0 everywhere, is equivalent to a ST theory, since in this case, the held equations solutions are 
the same of the usual ST (provided that the boundary conditions and the functions u (<h) and A (<h) 
are the same). Thence, the solutions of ST theory are particular solutions of the corresponding 
EST theory. However, ([ 5 ]) differs from the usual scalar equation by the scalar held gradient present 
as an overall factor. As a consequence, independently on the w (<f>) and A (#) functions, ([5]) also 


2 Similarly, for electromagnetism within an external spacetime geometry (that can be Minkowski in particular), 
one can use the action S[A] — — \ f F a bF ab ^/—g d 4 x to derive the field equations. These are derived by setting 
5S/8A a = 0. Note that the equation 5S/5g ab = 0 is wrong in general in this context. 

This is obtained by varying |l| with respect to <E>, where is not an external field in the usual ST framework. 


3 
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admits solutions such that 

d b $ = 0. (7) 

Reinserting in ([2]), it immediately results that the EST field equation ([2]) becomes identical to GR 
field equation, with the Newton’s constant G numerically given by 1/$. Thence, GR solutions also 
generate EST solutions with the same stress tensor, the EST metric being the GR’s one and the 
scalar field being $ = 1/G, where G is Newton’s constant. 

To be complete, it could be worth identifying the cases where EST gravity reduces to the usual 
ST case. The difference between the two theories just results from the presence of 9<3? as an overall 
factor in ([5]). Thence, the reduction to usual ST occurs if, and only if, the solutions of <9<h = 0 also 
solve eq. ©>• This conditions reads (considering finite u and ui') 

* 3 si(i) +4 " r = 0 - < 8) 

It is interesting to point out that, if no cosmological term enters the theory (i.e., if A = 0), then 
EST and ST are equivalent if T = 0. Hence, in particular, vacuum solutions of EST gravity 
with A(<f>) = 0 are always ST gravity solutions. For instance, in this case, the propagation of 
gravitational waves in vacuum obeys the same equations in both theories. 

Another issue to be stressed is (and this is the main point of this paper) is that the set of 
EST solutions is larger than the set of the corresponding ST or GR solutions. Indeed, one can 
also derive EST solutions by matching parts of usual ST solutions to parts of GR solutions, if the 
relevant conditions are satisfied. These conditions ensure that the EST field equations are valid at 
all the spacetime points (this will be exemplified in the next section). Reciprocally, it results from 
eq. ([5]) that the general EST solution can only be built from joined ST and GR solutions. Let us 
remark that, in an EST solution that contains two disjoint spacetime regions where = 0, the 
solutions in these regions satisfy the GR equation, but in general with different numerical values 
for the Newtonian constant. 


III. A COSMOLOGICAL SOLUTION 


Let us consider in this section a specific EST action characterized by A = 0 and a constant a 
The usual ST version of this case is the original BD theory [26] . Consider a flat dust filled RW 
Universe of metric 


ds 2 = — 


dt 2 + a (t) 2 5ijdx l dxG 


(9) 
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The EST equations to be satisfied by the scale factor a(t), the scalar field <1> (t) and the dust 
density e ( t ) are © and Q. The non-trivial Einstein’s equation components can be written as 

./2 


, a „ uj " a'. 

3<h^- = 8-rre + — —-3—T , 

a 2 2 $ a 


( 10 ) 


J! 


J2 


a" 4>" u ; T' 2 a' <!>' a" 

2— + — = -w^rw-2-— --5-. 


a 2 <h 2 a a* 

Where the prime represents derivation with respect to time. The stress tensor conservation yields 


(ea 3 )' = 0. 


(ID 


For the present case, eq. ([5]) reads 


(a 3 $')' - 


8irea 3 


= 0. 


( 12 ) 


2w + 3_ 

The above equation is implicitly contained in the system (10) and We stress that, within 

the usual BD framework, the scalar field equation ([6]) should be considered, and this equation is 
not compatible with = 0. Both the GR and BD solutions of the flat dust filled RW Universe 
are known, and are solutions of (10) and ( |11[ ), as it can easily be checked. After reminding these 
solutions in the first subsection and defining the useful related notations, we built in the second 
subsection a non-trivial EST RW solution got joining a part of the GR solution to a part of the 
BD one. 


A. General Relativity and Brans-Dicke like solutions 

The GR solution reads acR{t) = A (t — to) 2 ^ 3 , where A and to ar e integration constants, to 
being the GR big bang time. From the EST perspective, this solution can be achieved if the field 
is set to be a constant. Thence, we write the corresponding EST solution in the form 

a GR (t) = A(t- t 0 ) 2/3 , (13) 

$GR (t) = C, 

where C is a constant. 

The BD solution was found in m- It can be equivalently rewritten as 

a B D (t) = B (t - t+) q+ (t - t-) q ~ , (14) 

$bd (t) = 4> 0 (t — t + ) p+ (t — t-) p ~ 
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with 


1 + St/1 + 5 W 


q± = 


4 + 3w 


(15) 


P± = 


1 ± 3s^/l + gw 


4 + 3w 

In the above, B, <J>o, t- and t+ are (independent) integration constants, and s is a sign. Remark 
that we can decide that t+ > t- without loss of generality (permuting t + and f_ is equivalent to 
change s in — s). We make this choice in the following, with the consequence that t+ is the BD 
singular time for observers living in the t > t+ phas<0 Let us also remark that 

2 


a BD^BD ~ ^0 B 3 

while the scalar equation dfjl) reads 


4 + 3w 


t - ( p+t- +p-t +) 




87reGT 


3 + 2w 

These two equations yield, reminding that ea 3 is constant thanks to the energy conservation 

87 rea 3 24> n B 3 


(16) 


(17) 


(18) 


3 + 2w 4 + 3w 

For simplicity, let us just consider scenarios with e > 0 (that is always the case considering ordinary 
matter) and <J>o > 0 (meaning that the local effective gravitational constant G e // = 4+2iJ 1 


3+2o; $5£) 


IS 


positive). The solutions (14) that fulfil these requirements exist if and only if uj satifies 


4 + 3w > 0. 


(19) 


Let us point out that the limit for uj —> 00 of the BD solution (14) reads, if t+ and t_ are 


considered constant in the limit process (which makes sense since there is no constraint on t + and 


t- for (14) to be a BD solution), 

&BD lim (t) = B(t- t+) 1/3 (t - t_) 1/3 


( 20 ) 


4> 


BD lim 


it) = $ 


0- 


This does not correspond to (13) unless requiring t+ = t _As pointed in [28J this is due to the fact 

that, generically, the limit of BD dust cosmology for w —> 00 is not just dust filled GR cosmology, 


For observers living in the phase t > t+, this singular time corresponds to a big bang time whatever s in the case 
uj > 0, but only if s = — 1 in the case w < 0. 


4 
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but GR cosmology with a matter sector containing a massless scalar field besides the original dust 
field. The presence of such a residual massless scalar is generic in such a limit process: it results 
from the fact that the BD action term w$ _1 (5$) 2 has no trivial limit when oj —> oo, the scalar 


gradient (9<h going then to zero [ 251 [29] . Accordingly, the residual scalar field entering the solution 
(20) is zero in the case t+ = i_, and in this case onlyj^] 


B. A nontrivial External Scalar-Tensor solution 


We develop here a nontrivial cosmological EST solution with u > —4/3. This solution is 
composed by two phases, one which is part of a GR solution (at earlier times, eq. ©) , and the 
other a ST one (at late times, eq. (fl4|) ). The instant that joins both solutions is labelled t m . 
It should be stressed that in order for this junction to be an EST solution, eqs. (10) and © 
need to be satisfied at any instant, including the instant t rn . The whole problem then depends on 
the eight parameters (A,C,to,t m , B,$>o,t + ,t-). Since both a" and <h" explicitly enter the system 
(10) and ©> and since the matter distribution is regular, the continuity of the scale factor, the 
scalar field and their first derivatives is required. This constrains four of these quantities to be 
dependent of the four other ones. The scenario being a BD-like expansion following an earlier 
GR-expansion, it is natural to express the four BD parameters ( B , <l>o,t+,t—) as a function of the 
three GR parameters (A,C,to) and the matching time t m , that is left arbitrary. Remark that the 
invariance of the problem with respect to time translation allows to chose arbitrarily one of the 
instants entering the problem. It is natural to adopt the global solution singularity as the origin 
of time, i.e. to set to = 0. The matching time then satisfies t m > 0. It should also satisfy t m > t+, 
otherwise there is no matching and no place for the GR phase. 

Before proceeding, we remark that the proposed solution union may sound as a matching 
problem as discussed in mm, in particular the problem of matching a Schwarszchild spacetime 
region with an RW one [32/ 23], that results in the so called Einstein-Straus vacuole solution 
[25] , (For matching conditions in ST theories, see for instance [35].) It is worth pointing out the 
differences to the current problem. The metric and coordinates as written in Q represent the 
solution in the whole spacetime, and every solution of the system (10) and © is an EST solution 
in the whole spacetime. It is necessary and sufficient that the equations © and © are satisfied, 
in particular this implies that one should just verify if these equations are satisfied for t < t m , 

5 This point is not just formal, it is physically meaningful. For instance, the age of a dust filled flat RW Universe 
within GR reads Tgr = 2/ (3Ho). On the other hand, the age of a dust filled flat RW Universe within BD, resulting 
from (141 and lj — > oo differs from GR, since all the values in the interval ] 1/(3 JTq), 2/(3Hq)\ can be achieved 
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t > t m and at t = t m . The proposed solution satisfies this criteria. Within the present context, for 


t < t m the EST solution comes from eq. (13), and for t > t m from eq. (14). Hence, the related 
dust densities are given by 


6 vr egr (t) = 


( 21 ) 


4:TT€BD {t) = 


3 + 2cj 


4>o 


4 + 3 u(t- t + f q+ (t - t-f q - ' 

The continuities of the scalar and its (logarithmic) derivative, of the scale factor and its (logarith¬ 
mic) derivative, give respectively 


4> 0 ( t m -t + ) p+ = c, 


( 22 ) 


p+ +^- = o, 


tm. t -l tm t— 


(23) 


B (t m — t + ) q+ (t m ~ f-) 9 =AtH 3 , 


(24) 


Q+ + Q- 


tm -t+ tm t— 3 t r 


Since p + + p- = 2/ (4 + 3w), one gets from (23) 

4 + 3a; 

tm t-\- = 


P+ (t+ -t-). 


(25) 


(26) 


Since t m > t + and 4 + 3a; > 0, this condition implies that p+ < 0 since we imposed t + > t_. 
Thence the sign s is necessarily —1, and the q± and p± exponents read 


q± = 


P± = 


1 + a; ± \ /1 + ^a; 


4 + 3a; 


(27) 


1 + 3 \ /1 + oO; 


4 + 3a; 

and have the following signs: > 0 if and only if u > 0, q + > 0, p_ > 0 and p + < 0. 


The quantities t + and t_ can be derived from (23), (25) and (27). One gets 

1 

— / =tr 


1 ±3 a /1 + fa; 


(28) 
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From the above, and since u > —4/3, one sees that t- < 0. One also sees from (28) that t + < t m , 


meaning that the matching really occurs without any further condition on t m . The quantities B 


and 4>o follow from (22) and (24). One gets 



(9 + 6w)“^ , 


(29) 


(9 + 6cj) 4 +3“ t, 


4+3cj 


The conditions above guarantee the continuity of a, 4*, a' and at any time, including t m . These 
conditions are not sufficient to guarantee the existence of a" and 4>" at t m . Since these quantities 


appear in the EST eq. (10), it may seem that the latter terms must exist in order for deriving an 


EST solution. Nevertheless, one can achieve an EST solution even if those second derivatives do 
not exist, since they only appear in the particular combination 2a"/a + $"/<&. Thence, only the 
continuity of this quantity, or equivalently of (<ha 2 ) , is required for the solution to make sense. 
We carefully check these points in Appendix [Bj Also, one can easily check from (21) that the dust 
density e is continuous at t m . 

With the above, it was possible to derive a nontrivial EST solution. The imposed conditions 
constrained the model parameters, nonetheless the value of t m is still arbitrary. This means that 
there is no way to derive the transition instant t m from the matching conditions. However, if this 
transition instant happened in the past of our universe, then a physical imprint of such transition 
should be left in the cosmological data. That is, from such imprint one can in principle fix the 
transition instant t m . Hence, t m is part of the boundary conditions of the problem. 


The case of a perfect fluid having a linear barotropic equation of state is also tractable and 
leads to very similar results, as shown in Appendix [Cj 


To illustrate these results, let us consider the case uj = 0, and normalise the relevant quantities 


setting A = t m = 1 and C = 2. The BD big bang time is then t+ = 1/4, and the solution (13) and 


(14) then reads, using (27), (28) and (29) 


aGR 


(t) = t 2/3 and Tor (t) = 2, 


(30) 


2/1 2 t + - 

aBD (i ) = 7ir 4 and ® BD ® = vs /—t 


The transition from the GR solution to the ST solution is shown in Figure [T] 
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t t 


(a) Evolution of a and 1? from GR to ST. 


(b) Second derivatives and continuity. 


FIG. 1: The time evolution of the scale factor a and the scalar field in the left plot, and the 
evolution of their second derivatives in the right plot. The transition happens at the instant 
t m = 1, and a was normalized such that a(t m ) = 1. 


This illustrates that, for any t m , a solution can be derived by joining the (t < t m ) GR-like 
solution to the (t > t m ) BD-like solution. It should be noted that the converse is also possible, but 
it can only happen at a time where the scalar field experiments an extremum. In this case, the 
cosmological history begins by an era where a scalar field is acting, resulting in an expansion that 
occurs faster than the GR case. At t = t m , the scalar held is frozen, and the subsequent universe’s 
expansion is governed by GR. Such a possibility may be worth to be considered in the context of 
inflationary models. 

Considering the solution with (a, $) ( t. < t m ) = (a, <F ) Gi 2 an d (°> $) (t > t m ) = and 

in accordance with the analytical result presented in Appendix [Bj 


using the same numerical values as in Figure la, the continuity of (<ha 2 ) is shown in Figure 


lb 


IV. DISCUSSION 

We showed in this paper that giving an external status to the scalar held entering a ST theory 
increases the set of solutions of the theory. Indeed, both GR and BD solutions are EST solutions, 
and it is also possible to match pieces of GR solutions to BD ones in such a way to generate a 
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new class of solutions that neither belong to GR nor to ST, but satisfies the EST equations. This 


is explicitly illustrated by the cosmological example considered in Section III We stress that this 
new solution is not possible to achieve from pure BD gravity since the transition from non-GR 
to GR-like solutions, within BD, requires u> —>• oo, while for an EST theory the latter can be 
achieved for arbitrary a; values. 

While GR and ST RW cosmologies constitute well-posed Cauchy problems, this is false for EST. 
Indeed, if an EST cosmology is at a GR-like phase at a certain instant, it is impossible to predict 
when the change to a BD-like phase will happen. This means that it is generally necessary to 
provide more data than the standard data provided for a Cauchy problem, the transition time t m 
being a part of the boundary conditions to be required. 

It is nontrivial that ST gravity with an external scalar field, which we called EST gravity, has 
solutions that are neither GR nor pure ST, as we have shown here. ST gravities commonly are 
capable of generating GR solutions under certain limits, while EST gravity has an additional route 
capable of doing the transition between GR and non-GR local solutions. In this scenario, certain 
spacetime regions behave as GR solutions, others as ST solutions, while the complete spacetime 
satisfies the EST equations that are given by eqs. ©> and ([3]), or equivalently by eqs. ([2]) and 0. 
This kind of behaviour is akin to screening mechanisms of gravity [5145] , but the mechanism of EST 
gravity is of nondynamical nature, depending on the boundary conditions. We have here explored 


some cosmological configurations in which nontrivial EST solutions can be derived (Section III 


Appendix [0, and we expect to explore non homogeneous and more realistic configurations in a 
future publication. 

In this work, we have specially explored the use of an external scalar as an extension to Brans- 
Dicke-like actions 0 . In Sec. |III| explicit results were developed for the particular case of constant 
uj and no potential (A = 0). The extension of this analysis towards nonzero potential is surely 
important. For example, it would allow to apply the EST proposal to f(R ) theories in such a way 
that in some spacetime regions either GR or certain f(R ) theory would be valid. That is due to 
the fact that /(/?.) theories can be expressed by scalar-tensor theories with a non-null potential, 
and this is valid for diverse f(R) approaches, from the standard one, to Palatini, Hybrid gravity 
and others [391431] . Thus this EST application would yield an extension of f(R) gravity such that 
neither GR and nor standard f(R) gravity would be valid in the complete spacetime, but it would 
be possible to subdivide the spacetime in regions in which either one of the two would be valid. By 
using the usual correspondence between standard f(R) gravity and ST theories, it is not hard to 
see that this EST application to f(R) theories yields an f(R) extension in which all the possible 
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metric variations must be considered at the action level, except for those that change the Ricci 
scalar. Further developments on this EST application to /(/?.) theories should appear in a future 
work. 
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Appendix A: From diffeomorphism invariance to the stress tensor conservation 


Consider the diffeomorphism 


X ,a = x a +e ( x c ) 


(Al) 


If the metric’s components are considered as functions of the coordinates (rather than as functions 
of the point), the difference 5g ab = g' ab (x) — g a b(x) (that is not the components’ variation at a 


given point) induced by (Al) reads, at first order in £ [3TJ 


dgab = -gacdbZ c - gbcd a £ c - i c d c g a b- 


(A2) 


Following the same lines, any scalar held q i>, while its numerical value at any point is a gauge 


invariant, is varied as 5<j> = <j>' (x) — 4>(x) = —£ c <9 c </> by (Al) if cj) is regarded as a function of the 
coordinates. The diffeomorphism results in the following matter sector of ([Tj) variation 


5S m = 


d W 9^m) J.T, . d W gL m ) . \ , 4 


-<5'P + 


5g a b d x 


(A3) 


<5 4 ' 5g ab 

where the matter Lagrangian L m is defined by S m = J d 4 Xy/—gL m , and where the symbol jp 
represents the usual functional derivative of l w.r.t. the held F. If the part L m of the Lagrangian 
is a scalar, f d 4 Xy/—gL m is an invariant integral and 5S m is then necessarily zero whatever £. If 
none of the matter helds ’P is external, each T generates a Lagrange equation, and since all the 
\P occurrences in the total action ([I]) enters the matter sector, this Lagrange equation just reads 


S(V^gL 

in') _ 


<5\I> 


= 0. Thence, one gets, using the stress tensor dehnition T ab = 


ab _ 2 


\/~~g S 9a 


and (A2) 


J V^T ab (< g ac d b e + gM c + i c d c g ab ) d A x = 0. (A4) 

Since this occurs whatever £, one gets, after some usual algebra, the conservation equation ([3]). 
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Let us remark that ([3]) is not recovered if the scalar field 4? enters S m , for two reasons. The first 
one, relevant in the context of the present paper, is that T is an external field, and then does not 
generate a Lagrange equation. Second, even in the case were 4> would be non external, the resulting 
Lagrange equation = 0 would not reduce to = 0, since the matter sector does 

not contain all the scalar contributions to the action. 


Appendix B: On the continuity of ($c 


At t = t m , one can easily compute the left and right limits of the second derivatives, using (13) 


with to = 0, (14), (27), (28), (29) and the continuity relations (22), (23), (24), (25). For the scale 
factor, one gets 


2 A 

{t m ) = a GR (tm) = 4 /^: 

9 t m 


(Bl) 


Om) = a-BD (tm) = -■ 


4(3 + u)A 


9 (3 + 2w) 4 /3 ’ 

and they are different, meaning that a" is discontinuous at t m (thence that a" (t m ) does not exist). 
The same occurs for the scalar field 


(t m ) = &CR (tm) = 0, 


4 c 


<*>" ( 4 ) - &BD (tm) ~ 3 (3 + 2u) t 2 m 

Now, consider the quantity ($a?)". In the GR phase 

($a 2 ) GR = CA 2 t A/3 . 


Thence 


In the BD phase 


($a 2 )" (tm) = (^« 2 )gh (tm) = 


(<5>cl 2 ) bd = B 2 (f - t + ) 2q++p+ (t - t-) 2q - +v ~ 


Thence, using (25), 


(<ha 2 ) 


BD 


(4>a 2 


(tm) — „, 2 


16 2 q + +p + 2 (/_+p_ 


I BD 


9t m (t m t_||_) (t m t_) 


and one gets 


($a 2 )" (t+) = ($a 2 )” BD {t m ) = = (4>a 2 )" (t m ) . 

9/ rn 


(B2) 


(B3) 


(B4) 


(B5) 


(B6) 


(B7) 


Since (4>o 2 ) / (t m ) is well-defined, (Qa 2 )" is then continuous at t m , with the value just calculated. 
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Appendix C: The general linear barotropic case 


In Section |III B it is presented a cosmological model whose matter content is a simple dust 
fluid. This simple kind of matter was used since the main purpose was to show the existence of 
nontrivial EST solutions. In this appendix we show that the existence of nontrivial EST solutions 
can also be found in systems with non negligible pressure, in particular we consider a barotropic 
fluid. It turns out that the case of a perfect fluid with a barotropic equation of state having the 
form P = ne is also tractable, leading to similar computations. Indeed, the BD solution reads in 
this case |27J (just considering theories with u > —3/2) 

a>BD (t) = B(t- t + ) Q+ {t- t-) Q ~ , 


®bd (t) = $0 (t — t + ) + {t — t_) 


£bd ( t ) clbd (t) 


3+3n 


= M 


(Cl) 

(C2) 

(C3) 


with 


Q± = 

1 H - (1 — 71 ) UJ zb yj 1 + 

(C4) 

4 — 6n + 3 (1 — n) 2 uj 

P± = 

1 — 3n =F 3 (1 — n) yjl + 

(C5) 

4 — 6n + 3 (1 — n) 2 oj ’ 


where t is a conformal time defined by 


dt = a 3n dt. 


(C6) 


The constants B, $ 0 ) t + and M are the five integration constants. The GR solution acR (t) oc 
(t — t 0 ) 2/(3+3n) can be written, using the same conformal time, and choosing the origin of time 
such that to = 0, 

a GR (t) = At i ^ 3 - 3n \ (C7) 

$gr (t) = C, (C8) 

e GR (t) OLGR (t) 3+3ri = M. (C9) 

The continuity of the scale factor and the scalar and of their derivatives with respect to t result in 

dt d 
dt dt 


the continuity of their derivatives with respect to the conformal time t, since = a3n J;- 


These continuities at t m lead to a system that can be put on a form very close to (22)-(25), 

V P + 


$o {tm - t+) + {tm - t~) = C, 


(CIO) 
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P- 


tm, 


+ 


P_ 


i'rn. 


= 0 , 


(Cll) 


B (t m -t + f~ n)Q+ (t m -t.f- n)Q - = AVi\ 


(C12) 


(1 — n)Q+ (l — n)Q-_ 2 


(C13) 


The solution for the BD constants is then the same as for the dust case, just replacing the exponents 
p± and q± by P± and (1 — n) Q± respectively. 


Appendix D: A mechanical example with an external coordinate 


This appendix aims to clarify the physical meaning of external fields by providing an additional 
example within a very simple context, that of a planar double pendulum. 

Before considering the case with external coordinates, consider the usual Lagrangian for the 
double pendulum subject to a uniform gravitational acceleration whose norm is denoted by g ( 36| . 

L(6i,0i ) = ^d\l\ (mi + m 2 ) + ^0fZ|m 2 + 9i9 2 hl 2 m 2 cos (9\ - 9 2 ) - V(9i), (Dl) 

with 


V ( 9i ) = -ghm 2 cos (6b) - gli (mi + m 2 ) cos (9 1 ), 


(D2) 


where 9i is used as a shorthand notation for 9i,9 2 , the particle masses are denoted by m\ and m 2 , 
mi is supported by a wire of negligible mass of length l\, m\ and m 2 are connected by another 
wire of length l 2 , 9\ and 9 2 are the position angles of each of the particles in regard to a straight 
line parallel to the gravitational acceleration vector g. 

The equations of motion read 


m 2 \ l\ •• 


mi ) l 2 


m 2 


^2a 2 


9 


1 H- — 9 \ H- 9 2 cos (6*i — 9 2 ) -\ - 9 % sin { 9 \ — 9 2 ) + — 1 H-sin (6>i) = 0, (D3) 


mi 


m i 


h 


m 2 


m i 




+ 6i cos (6*i — 9 2 ) — 9\ sin ( 9\ - 9 2 ) + 


j- sin (0 2 ) = 0. 
n 


(D4) 
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Considering that rn \ is sufficiently large, the above equations can be written as 

hOi =-gsm(6i), (D5) 

y-0 2 = - ^ cos (Q l _ q 2 ) + Qj s i n [Q x - 0 2 ). (D6) 

£1 £1 

In this limit, the second pendulum has no effect on the first one. Alternatively, one can also consider 
the regime rn 1 m 2 right in the Lagrangian (or the action), and the resulting Lagrangian reads 

L(6i, 9i) = m 2 | + 9\9 2 l\h COS (9± - 9 2 ) + gl 2 cos(02)| + 

+ (mi + m 2 ) (gh cos (0i) + 

= m,2L 2 {9i,di) + (mi + m 2 )Li(0i, 0i) 


m 2 L 2 (92,92,t) +miL l (9i,9i). 


(D7) 


The last approximation is constituted by two parts, the first is a trivial one in which it was used 
that mi + m 2 ~ mi. The second part refers to the replacement L 2 (9i,9i) —> L 2 (9 2 ,9 2 ,t), in which 
L 2 becomes a function of the second particle only, and the occurrence of 9\ inside L 2 is seen just 
as function of the time variable t, not a generalised coordinate with the same status of 0 2 . This 
type of decoupling on the dynamics of the particles 1 and 2 is expected since, by requiring that 
mi is sufficiently large, the extremum of the total action f Ldt should be achieved only when 
both m 2 f L 2 dt and mi / L\dt are independently extremized. Hence, a consistent approach is the 
following: from the Euler-Lagrange equations of Li one derives #i(t); and with this result, the 9 1 
solution is inserted into L 2 , which is then used to derive the 9 2 solution. 

Another way of justifying why it is correct and necessary to use L 2 (9i,9i) —» L 2 (9 2 ,9 2 ,t), 
when the regime mi 3> m 2 is considered, comes from the fact that the Euler-Lagrange equations 


of L\(9i,9\) leads to the eq. (D5), while the Euler-Lagrange equation of L 2 (9 2 ,9 2 ,t) leads to eq. 


(D6). Moreover, if one considers L 2 as a function of 9\ and derive its corresponding Euler-Lagrange 


equation, the derived equation is not compatible with eq. (D6). Hence it is wrong to deal with 9\ 


in L 2 as if it were a standard coordinate: 9 1 is an external coordinate in L 2 (9 2 ,9 2 ,t). 


In the following a generalisation of the picture above presented will be considered. Eventually, the 
dynamics of the external coordinate may not be well known, that is, the Lagrangian Li may be 


more complex. Since eq. (D6) does not depend on Li, one can use L 2 in a more general situation 
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in which L\ is unknown. In this setting, it will not be possible in general to derive the full solution 
of both 6 \ and 82 , unless additional information is provided. Nonetheless we stress that it is correct 


to use eq. (D6) even though L 1 , and consequently is unknown. 

Consider the Lagrangian 

L(0i,6i) = L(8i, Oi) + kL 3 { 8 1 , Mi, ■■•)> 


(D8) 


where L is given by eq. (Dl), k is a dimensionless quantity, and L 3 is any function that depends 
only on 8 \ and its derivatives (it may or may not depend on higher derivatives). Here it is not 
assumed that m\ m 2 . Nonetheless, for k sufficiently large the same argument used in the case 
in which mi is sufficiently large can be used. Therefore the extremum of the action f Ldt will be 
achieved only when both f Ldt and f L 3 dt are extremized, thus, in the regime k 1, one can 
write 


L(8i,8 t )^L(02,82,t) + kL 3 (e (D9) 

where the explicit time dependence of L comes from the solution of 9i(t), which in turn can be 
directly derived from the equations of motion of / L 3 dt. It should be noted that even if the 
Lagrangian L 3 is unknown, one can use the equation of motion of L 2 with respect to 82 , while 8 \ 
is seen as an external coordinate. In this case, 8 \ (t) will be an arbitrary function. 


Appendix E: A first test of an EST model with SN data 

The main intention of the present paper is to show the existence of nontrivial solutions within 
the EST approach. It has been shown in other sections that the matching between GR and BD 
solutions, in the context of the EST approach, is possible. Although, for the moment, the resulting 
cosmological model is too simplified to constitute a realistic scenario, it is already possible to test 
the configuration here obtained with SN data. This is done here both to show that it is possible 
and as a first step towards testing EST models. To this end, we use the binned JLA sample m- 
The advantage of such sample is related to the binned process which lead to less important effects 
due to different calibration methods. 

We match the flat CDM model of GR with the flat BD solutions used above. Of course, 
the observational test disfavour strongly the CDM with respect to the ACDM model. We use the 
expressions derived above, and consider three free parameters, the matching time t m , the parameter 
oj and the Hubble parameter h. As usual, the fit is done looking for the values of these parameters 
that lead to the global minimum of \' 2 . 
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The results indicate that the global maximum of probability (or global minimum of y 2 ) occurs 
for t m = to, where to is the present time, implying that the data favours no phase transition and the 
universe is still in the GR phase. The derived value of the minimum \' 2 is not small in comparison 
to the standard cosmological model ACDM, but this is expected since this model includes no dark 
energy. The resulting fit is just as good as that of a pure CDM model. Likewise pure CDM, the 
minimum x 2 , for this set of data, is x 2 ~ 350. There is a local maximum for t m = 0 corresponding 
to a pure BD phase, for the latter, the corresponding value is % 2 ~ 353. For comparison purposes, 
the minimum x 2 for ACDM is 29. The la and 2 a confidence regions are displayed in fig. ([ 2 ]). As 
expected, the higher is u, more similar the dynamics of the BD and the GR phases become, and 
hence it becomes more probable for the transition to happen in an instant before to- 

A complete statistical study would imply to consider other tests and to introduce a mechanism 
to accelerate the expansion of the universe. This has in part been done in Ref. [38] in the context 
of the usual BD theory by including the possibility of w < —4/3, what has been excluded in our 
analysis. In that reference, it has been found, comparing the pure BD flat matter dominated phase 
with the ACDM, that the BD model is slightly favoured for uj < —4/3. 



OJ 


FIG. 2: The la and 2a confidence regions in the ( uj,t m ) space. The SN fit was done with three 
free parameters, the latter two and the Hubble parameter ( h ). For the result above, h was 

marginalised. 
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